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This  report  documents  one  of  the  perturbing  terms  appearing  in  the  force 
equations  for  satellite  motion  in  the  TERRA  system  of  satellite  geodesy  computer 
programs.  Equations  are  given  for  the  lunar  air  tide,  via  the  gravitational  action 
which  the  tidal  bulge  exerts  on  the  satellite  orbit.  An  earlier  report  gave  the 
equations  for  the  air  tides  caused  by  solar  gravitational  and  thermal  effects. 

The  work  was  done  in  the  Astronautics  and  Geodesy  Division  in  support  of 
astronautics  computer  program  development. 
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INTRODUCTION 
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The  atmospheric  tides  of  the  earth  are  caused  by  a combination  of  the 

gravitational  pull  of  the  moon  and  the  sun’s  thermal  action.  To  establish  a 
mathematical  model  of  the  tidal  effects  on  satellite  motion,  it  is  both  convenient 
and  customary1  to  represent  the  tidal  mass  redistribution  by  a surface  density  layer 
and  to  subsequently  obtain  the  potential,  exterior  to  the  earth,  by  integration. 

One  of  the  authors  has  recently  derived  the  solar  component  of  the 
atmospheric  tide’  from  a surface  pressure  pattern  described  by  Haurwitz.1  This  was 
done  to  produce  a perturbing  term,  corresponding  to  the  solar  air  tide,  for  the 

equations  of  motion  which  govern  the  orbits  of  the  earth  satellites  in  the  TERRA1 
system  of  computer  programs  for  satellite  geodesy.  1 or  the  same  purpose,  the 

following  text  introduces  a model  which  accounts  for  the  lunar  tide  in  the 

terrestrial  atmosphere.  Specifically,  the  new  tide  model  contains  a potential  function 
for  the  force  by  which  the  atmospheric  tide  bulge  acts  on  the  satellite  orbits.  The 

tidal  bulge  is  assumed  to  result  from  the  fact  that  the  earth  rotates  within  the  field 

of  lunar  mass  attraction,  the  latter  being  inhomogeneous  across  the  terrestrial  globe. 

Only  the  main  term  of  the  semidiurnal  tide  is  considered.  The  perturbing 
acceleration  associated  with  the  tide  potential  is  also  specified  for  use  with  TERRA 
and  with  similar  computer  programs  which  involve  force  equations  from  which  the 
orbits  of  geodetic  satellites  are  being  calculated  by  direct,  numerical  integration. 

i] 

The  text  below  augments  Reference  2.  Readers  ought  to  have  peruseil  that 

reference  first,  especially  the  introduction.  Much  has  been  included  there  to  facilitate 
the  understanding  of  the  more  general  aspects  of  the  matter,  which  we  shall  not 

repeat  here 


THE  TIDE  POTENTIAL* 


A suitable  description  of  the  variation  in  air  pressure  on  the  ground,  caused  by 
the  semidiurnal  lunar  tide,  was  found  in  Reference  4 (pages  70,  73.  93,  and  94): 


L,  = 55.2jubar  cos'  (I  sin  (2r  + 75°) 


(I) 


* Additional  comments  on  the  tide  potenti.il  plus  the  necessary  details  of  its  derivation  appear  in  Appendix  A. 


1 


In  this  equation,  L,  is  the  surface  pressure  in  inicrobars,  0 is  the  geocentric 
latitude,  and  t is  the  local,  lunar  mean  time.  Considering  that  the  local  surface 
pressure  results  from  the  weight  of  the  atmospheric  column,  and  also  taking  into 
account  that  geodetic  satellites  generally  operate  far  above  that  altitude  below  which 
most  of  the  atmosphere  lies,  the  surface  pressure  was  converted  into  a gravitating 
mass  layer,  the  action  of  which  on  the  satellite  would  resemble  that  of  the  actual 
tidal  air  bulge: 


= A,  cos3  0 cos  2a 

(2) 

= 0.0564  gr/cnr 

(3) 

2\  + [2(t*  - v)  - 1 5°  1 

(4) 

In  these  equations,  X is  longitude  bast,  t*  is  Universal  time,  and  v is  the  difference 
between  the  mean  mean  longitudes  of  moon  and  son. 

Subsequent  application  of  the  Poisson  integral  produced  the  lunar  tide  potential, 
exterior  to  the  main  portion  of  the  atmosphere. 


R3 
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P;(sin  0)  cos  2 a - b — 
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Cl  is  the  gravitational  constant,  R is  a suitably  defined  mean  radius  of  the  earth, 
and  r is  the  geocentric  distance  of  the  satellite.  The  P -Jl x ) are  associated  Legendre 
functions. 

It  should  be  noted  that  L,  above  is  a yearly,  worldwide  average.  Also,  all 
equations  appearing  here  and  in  the  text  below  are  presented  in  a form  which 


renders  them  independent  of  the  choice 
compatible  units  must,  however,  be  chosen 
into  a computer  algorithm. 


of  any  particular  units.  A system  of 
prior  to  translating  the  present  analysis 


THE-  MEAN  LONGITUDES  OE  SUN  AND  MOON* 


Explicitly  and  implicitly,  the  local,  lunar  mean  time  t is  the  measure  of  time 
in  the  equations  of  the  previous  chapter.  When  converted  into  Universal  time,  it 
gave  rise  to  the  presence  of  the  quantity  v in  the  combined  azimuth/time  arguments 
of  surface  density  and  potential. 

v is  the  difference  between  the  mean  mean  longitude  of  the  moon,  s,  and  the 
mean  mean  longitude  of  the  sun,  h: 


v = s - h 


(8) 


To  be  precise,  s is  the  angle  front  the  mean  equinox  of  date,  along  the  mean 
ecliptic  of  date,  to  the  mean  ascending  node  of  the  moon’s  orbit,  plus  the  angle 
from  that  latter  point,  along  the  moon’s  mean  orbit,  to  the  mean  perigee,  plus  the 
instantaneous  value  of  the  mean  mean  anomaly  corresponding  to  the  position  of  the 
moon.  Correspondingly,  h is  the  angle  from  the  mean  equinox  of  date,  along  the 
mean  ecliptic  of  date,  to  the  mean  solar  perigee,  plus  the  instantaneous  value  of  the 
mean  mean  anomaly  corresponding  to  the  position  of  the  sun. 

The  terms,  s and  It.  are  connected  to  Universal  time  t*  by  the  following 
algorithm.  Let  J be  the  number  of  the  calendar  year  ( 1 ‘>75  or  l‘)7(>  or  W77),  M 
be  the  number  of  the  individual  month  in  the  calendar  year,  and  1)  be  the  number 
of  the  particular  day  in  the  month.  Also,  observe 

II  n = m 

if  (9) 

0 n ¥=  m 


See  also  Appendix  11. 


Then 


N,  = 0 + 3 15(M,2)  + 595(M,3)  + 906(M,4)  + 1205(M,5)+  1515(M,6) 

+ 1 8 1 fi(  M,7)  + 21 26(M,8)  + 2436(M,9)  + 2735(M,10) 

+ 3045(M,1 1 ) + 3346(M,1 2)  (10) 


N,  = 1>  + 315(M,2)  + 605(M,3)  + 91S(M,4)  + 1215(M,5)+  1526(M.6) 

+ 1825(M,7)  + 2138(M,8)  + 2445(M,9)  + 2745(M,10) 

+ 3055(M,  1 1 ) + 3355(M , 1 2)  (11) 

N = N,  5(J,1975)  + (365  + N,  )6( J, 1 976)  + (731  + N,  )fi(J,1977)  (12) 

AT  = 5.28  X 10  4 +3.56  X 10  8 X N (13) 

t* 

d = 27392.5  + N + AT  + — (14) 

24 


(Into  this  equation,  t*  should  be  entered  in  decimal  fractions  of  hours.) 


d 

r = 

36525 


(15) 


s = 270.434358  + 48 1267.883 14 IT-  0.001 133T2  + 0.000002T1  ( lo) 


h = 279.69668  + 36000. 768930T  + 0.000303T2 


(17) 


(s  and  h will  result  in  decimal  fractions  of  degrees.) 

It  is  obvious  how  this  algorithm  may  be  generalized  to  time  spans  which 
extend  beyond  1977  or  to  altogether  different  time  spans.  Also,  one  may  desire  to 
input  the  number  of  the  day  in  the  calendar  year,  N(  or  N,,  directly.  That  is 
easily  accomplished  by  deleting  the  equations  for  Nj  and  N,.  Nf  and  N,  would 


4 


then  be  algorithm  input,  to  be  obtained  from  an  external  calendar.  Nj  would  then 
mean  the  number  of  the  calendar  day  in  the  year,  if  J = 1975  or  J = 1977.  N., 
would  be  the  number  of  the  calendar  day,  if  J = 1976. 


DISTURBING  ACCELERATION* 

From  the  potential  U (Equation  5),  the  disturbing  acceleration  T can  now 
readily  be  obtained.  A positive  gradient  convention  will  be  observed, 

T = + grad  U (18) 

Let  y ! , y.,,  y,  be  the  earth-fixed,  Cartesian  coordinates  of  the  TERRA  system. 
Also,  let  r,  X,  0 be  the  above  used,  earth-fixed,  spherical  coordinates  (geocentric 
distance,  longitude  East,  and  geocentric  latitude).  Further,  observe  the  geometrical 
relations. 


r = 

^y;*y;+y; 

(19) 

tan  X = 

V| 

(20) 

cos  X = 

yi 

(21) 

sin  X = 

y? 

(22) 

^ y;  + y; 

tan  0 = 

y< 

(23) 

>1  y]  + V*: 


*SlV  ;lln>  Appi'IKllN  ('. 
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L. 


cos  0 

= 

1 n“ : 2 

-^y-  + y; 

(24) 

r l 1 

sin  0 

= 

y3/r 

(25) 

Or 

y= 

9Vi 

Ll,  i = 1,2.3 

r 

(26) 

ax 

y2 

(27) 

9V| 

y? + V2 

ax 

. y> 

(28) 

ay: 

+ ■>  , 2 
yy  + y3 

ax 

0 

9y3 

(29) 

be 

- yi  y3 

(30) 

3y, 

r>Jy,  + y\ 

bO 

- y2  y,3 

(31) 

r'j 

>> 

rz> 

2 1"  2 2 

r\yl  + y2 

bO 

•^yy  + y^ 

(32) 

>» 

ro 

r2 

Also  observe  the  physical  relations. 


9U  3R^  ->r5 

T — = -a — t-  Pilsin  0)  cos  2a  + b z~  P: (sin  0)  cos  2a 
3r  r 2 r°  4 


(33) 


9U  2R3  , 2R5  , 

tt  - -a — r-  P2  (sin  0)  sin  2a  + b— Pr  (sin  0)  sin  2a 

3X  r 2 r5  4 


(34) 
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9 1.>  j>i\  i.ii\ 

- — = -a— 5-  sin  20  cos  2a  + b — 7- 
dO  r3  r5 


(7  sin3  0-4)  sin  20  cos  2a 


P;(x)  = 3(1  - x2) 

Pj(x)  = ^(1  - x2)(7x2  - I) 


x = sin  0 


_yi 

r 


sin  20 


cos  2a 


•> 


COS' 


a - sin2  a 


sin  2a  = 2 sin  a cos  a 


a = X + a* 


a*  = t*  - v - 7.5° 


v = s - h 


cos  a = cos  a*  cos  X - sin  a*  sin  X 


sin  a = sin  a*  cos  X + cos  a*  sin  X 


Then,  the  Cartesian  components  of  the  disturbing  acceleration  are 


auar_  au  9u  a o 

9r  9y , 3X  9y,  9 0 3y ( 


(46a) 


9U  9r  9U  3X  9U  90 

j = + + (46b) 

y2  9r  3y,  9X  3y,  9 0 9y, 


9U_9i_  9U  9X_  3U  9 0 

•3  ” 3r  9y3  9X  9y3  3 0 9y3 


(46c) 


To  be  suitable  for  use  in  the  satellite  orbit  integration,  T will  now  have  to 
be  rotated  into  the  inertial  coordinate  frame  in  terms  of  which  IT' R RA  executes 
the  orbit  integration.  The  necessary  algorithm  shall  not  be  discussed  here.  It  is 
available  for  implementation  on  the  computer  now  and  will  appear  in  a future 
report  which  will  document  the  complete  TFRRA  equations  of  motion,  strictly 
under  an  algorithm  aspect. 

Terms  like  “perturbing  force”  and  “perturbing”  or  “disturbing  potential”  are 
frequently  associated  with  General  Theory.  To  prevent  any  misunderstanding,  it 
should  be  remembered  that  T1RRA.  like  its  forerunners,  generates  the  necessary 
trajectories  by  directly  integrating  a set  of  differential  equations  which  relate  the 
orbital  acceleration  to  a variety  of  force  terms,  which  in  turn  are  rather  complicated 
functions  of  satellite  position,  velocity,  attitude,  position  relative  to  the  earth’s 
shadow,  etc.  This  approach  enables  TFRRA  to  readily  cope  with  force  models  which 
reflect  the  state  of  the  art  and  which  are  mathematically  complex  to  the  point 
where  they  might  exceed  the  capabilities  of  General  Theory.  That  is  why  we  were 
able  to  formulate  the  algorithm  in  the  simple-minded  shape  in  which  it  appears,  and 
this  is  also  the  reason  why  it  lacks  the  elegance  which  is  so  characteristic  ot  the 
analytical  approach. 


SUMMARY 
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The  above  text  presents  the  essential  features  of  only  one  from  amongst  four 
tidal  force  terms  which  the  T1KRA  equations  of  motion  contain.  The  analytical 
work  has  been  completed  for  all  these  terms.  Also,  the  necessary  computer 
algorithms  have  been  formulated.  But  only  two  perturbing  terms  have  thus  far  been 
described  formally.  This  includes  the  solar  air  tide  which  is  covered  by  Reference  2 
and  the  present  attempt  to  elaborate  that  component  of  the  perturbing  acceleration 
which  arises  because  of  the  presence  of  the  lunar  tide  in  the  atmosphere.  1 wo  more 
reports,  namely  one  on  the  tide  of  the  solid  earth  as  well  as  another  covering  the 
ocean  tide,  are  in  preparation. 

Tile  T1RRA  equations  of  motion  are  intended  to  serve  as  a basis  lor  numerical 
orbit  integration.  While  developing  the  individual  force  terms,  we  never  gave  any 
thought  to  General  Theory.  In  keeping  with  this  approach  and  for  other,  entirely 
pragmatic  reasons  as  well,  it  was  decided  to  make  the  tidal  terms  part  of  Tl'RRA 
without  first  conducting  a study  of  the  relative  magnitudes  of  the  various  force 
components  and  their  consequences  for  the  satellite  motion.  That  might  have 
involved  us  in  orbital  studies  based  on  General  Theory,  which  was  considered 
undesirable  in  view  of  both  the  funding  situation  and  the  time  frame  of  the 
T1  RRA  development.  Rather,  it  was  felt  that  we  already  had  reliable  indications 
from  experience  as  well  as  from  the  literature  as  to  which  effects  could  be  expected 
to  be  relevant  and  which  might  safely  be  omitted,  l or  that  reason  and  also  because 
it  appeared  to  be  more  in  keeping  with  the  rest  of  the  T!  RRA  effort,  we  decided 
that  all  the  selected  forces  should  appear  in  the  equations  of  motion  exactly  as 
they  resulted  from  the  respective  physical  models  and  that  their  existence  should  be 
finally  justified  by  means  of  numerical  studies  which  would  have  to  be  conducted 
anyway,  as  part  of  the  Tl'RRA  trials  on  the  computer. 
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APPENDIX  A 

NOTES  ON  THE  TIDE  POTENTIAL 


NOUS  ON  THE  TIDE  POTENTIAL 


Our  work  on  the  lunar  air  tide  is  based  on  a tidal  pressure  wave  (Equation  1) 
which  is  symmetric  to  the  equator.  Also,  time  and  latitude  are  altogether  decoupled 
from  each  other.  Time  is  associated  with  the  azimuth  coordinate  only.  That  is,  of 
course,  a simplification.  Being  caused  by  the  moon,  the  orbit  of  which  is  inclined 
against  the  equator,  the  actual  tide  may  be  expected  to  depend  on  the  lunar  phase. 
That  excludes  strict,  equatorial  symmetry.  But  the  earth’s  topography  appears  to 
complicate  the  tidal  pressure  wave.  Also,  data  are  still  scarce  over  much  of  the 
globe.  After  having  reviewed  the  recent  literature,  we  were  quite  content  with  having 
available  at  last  something  tangible,  namely,  the  tidal  pressure  averages  represented 
by  Equation  1. 

The  surface  pressure  wave  we  used  can  be  imagined  to  be  caused  by  a moon 
which  moves  in  the  celestial  equator  rather  than  in  the  actual  lunar  orbit. 
Figure  A-l  indicates  how  the  various  azimuth  coordinates  appearing  above  in  the 
main  part  of  the  report,  and  time,  are  related  to  each  other. 


Figure  A-1.  Relations  Amongst  the  Azimuth  Variables  and  Time 


A-l 


t is  local  mean  solar  time,  t is  local  mean  lunar  time.  Also,  consider  the 
familiar  relationship  between  Universal  time  and  t, 


t*  = t-  X 


Then, 


r = t-i>  = t*  + X-*»  (48) 

which  leads  directly  to  Equation  4 for  the  time/azimuth  argument  of  the  tidal  mass 
layer. 

From  Equation  1 for  the  surface  pressure  follows  the  corresponding  surface 
mass  density  o,  by  simply  replacing  the  unit  of  pressure  by  the  area  density  related 
to  the  atmospheric  column.  The  potential  aloft  follows  by  integration: 


= f C C,q2(Q 

3 3 if- pi 


where  P and  Q denote  the  field  point  and  the  source  point,  respectively,  r is  the 
position  vector  for  the  field  point,  p indicates  the  source  point.  When  used  in  the 
following,  the  field  point  coordinates  will  be  unprimed  while  the  source  point 
coordinates  will  be  primed,  y is  the  angle  between  r and  p.  The  area  element  ds  is 


ds  = R-  cos  O’dO'da' 


There  will  appear  below  several  expansions  involving  Legendre  polynomials.  These 
were  taken  from  Reference  5. 

To  actually  perform  the  integration,  the  inverse  of  |F  p|  is  first  expanded  in 
terms  of  Legendre  polynomials: 


|r-  p|  = +-\r*+p  - 2rp  cos  7 


cos  7 = sin  0 sin  O'  + cos  0 cos  O'  cos  (a  a') 


(52) 
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IF  pi 
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“ 3Tm  Pn(cos7) 

n = 0 r 


(53) 


Application  of  the  Addition  Theorem  for  the  Legendre  functions  (Reference  5, 
Chapter  VII,  Paragraph  (>)  results  then  in 


II 


A , ( i R 
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Rn>l 
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(55) 


I = P (sin  0)P  (sin  O')  + 2 - 

1 n n 

m I 


(n  - m)! 
(n  + in)! 


P"*  (sin  C )PJ"  ( sin  O')  cos  ni(o  o') 

(50) 


If  one  now  integrates  over  the  a/imuth  coordinate  first,  it  will  at  once  he  realized 
that  I is  a fairly  simple  expression: 
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K)n  = I cos  m(a  - o')  cos  2o 'do' 


(57) 


(58) 


Inspection  of  K will  show  that  only  K,  is  non-zero. 


Km  = K, 5(m.  2) 


( 5l>) 


K , = n cos  2o 


(<>()) 


I,  wc  shall  bo  able  to  ignore  for  reasons  indicated  on  Page  A-5  of  Reference  2. 
1 Inis,  there  remains 
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II 
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(n-  2)! 
(n  + 2)! 


1*“ (sin  0) 


P* (sin  O'k- os4  0\W' 


(61) 


I xecution  ot  the  latter  integral  and  search  for  the  non-zero  I shows  that  I,  is 
zero  while  I,  and  I4  result  in  expressions  which,  upon  entering"  into  liquation  54, 
lead  to  the  potential  represented  by  liquations  5,  6,  and  7. 
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Appendix  A indicated  how  the  perturbing  potential  was  derived  from  the  tidal 
pressure  wave.  It  did  not,  however,  clarity  a problem  hidden  within  the  azimuth 

coordinate,  namely  that  of  how  to  reference  the  moon's  position  to  that  of  the 

sun.  That  question  is  solved  by  the  algorithm  for  the  mean  longitudes  of  sun  and 

moon,  specified  above  in  the  main  part  of  the  report.  We  shall  now  discuss  some 

details  which  concern  the  computation  of  the  solar  and  lunar  mean  longitudes. 

In  essence,  the  algorithm  is  based  on  Iquations  1(>  and  l7.  which  represent  the 
mean  mean  longitudes  of  sun  and  moon  as  functions  ol  time,  lhese  equations  are 
from  Pages  and  107  of  Reference  b.  I hey  are  part  ol  the  body  ol  astronomical 
relationships  on  which  the  Nautical  Almanac  bases  its  fundamental  ephemerides  ol 

sun  and  moon.  The  mean  mean  longitudes  s and  h are  polynomials  ol  I.  1 is  time 
in  Julian  centuries.  These  polynomials  are  valid  tor  the  present  century  d as  defined 
by  I q nation  14  is  time  in  days;  llKK)  January  0.5  I I (I  phemeris  time)  is  the 

epoch.  N is  the  number  of  days  elapsed  since  the  beginning  ol  the  calendar  year 

1075.  AT  is  the  Reduction  of  lime  which  is  the  difference  between  I- phemeris  time 
and  Universal  time.  Our  1 quation  13  actually  is  only  an  approximation  to  that 

difference,  chosen  to  be  valid  for  any  time  interval  during  which  we  anticipate  that 
it  may  be  required. 

As  already  mentioned,  epoch  for  counting  time  is  the  beginning  ol  our  century 
or,  to  be  precise,  the  start  of  the  calendar  year  1 000.  Although  the  concept  of  the 
Julian  day  is  simple,  its  exact  connection  with  other  time  units  such  as  the  calendar 
day  may  likely  confuse  those  who  only  occasionally  encounter  that  particular  unit. 
We  add  thus  three  illustrations,  namely  Figures  15-1.  15-2,  and  15-3.  I he  additional 
Figure  15-4  illustrates  the  role  of  the  Reduction  of  Time  in  the  relationship  between 
I phemeris  time  and  Universal  time. 
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NOTES  CONCl  RNING  THE  GRADIENT 


The  perturbing  acceleration,  as  discussed  above  in  the  main  portion  of  the 
icport.  will  not  require  much  additional  elaboration.  But  we  probably  should 
mention  that  the  integrations  which  resulted  in  the  tide  potential  as  well  as  the 
gradient  formation  itself  were  done  by  both  authors,  independently  of  each  other, 
with  great  care.  Further  confidence  in  the  results  was  derived  from  a check  on  the 
sign  of  the  gradient.  This  was  to  assure  that,  for  a low  orbit,  the  radial  gradient 
component  associated  with  any  individual  potential  term  actually  does  point  down 
(negative  sign)  towards  the  surface  whenever  the  satellite  happens  to  be  positioned 
over  a region  of  the  surface  where  the  corresponding  atmospheric  density  term  is 
positive  and  thus  attracts  rather  than  repels  the  satellite. 


C-l 


DISTRIBUTION 


Defense  Mapping  Agency,  Headquarters 
Naval  Observatory  Building  56 
Washington,  DC  20360 

Defense  Mapping  Agency 
Aerospace  Center 
St.  Louis.  MO  63118 
ATTN:  Richard  L.  Ealum 

Defense  Mapping  Agency 
Hydrographic  Center 
Washington,  DC  20390 

Defense  Mapping  Agency 
Topographic  Center 
6500  Brooks  Lane 
Washington,  DC  20315 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  VA  22314 

Chief  of  Naval  Research 
Department  of  the  Navy 
Arlington,  VA  22217 

Department  of  the  Air  Force 
HQ  Space  & Missile  Test  Center  (ALSO 
Vandenberg  Air  Force  Base,  CA  93437 
ATTN:  Sven  C.  Hrnberg 

Naval  Research  Laboratory 
Washington,  DC  20360 
ATTN : Code  4130 

U.  S.  Navy  Astronautics  Croup 
Point  Mugu,  CA  93041 
ATTN : Code  AC, 33 

Code  SPMOO 
Mr.  Campbell 


Commanding  Officer 

Naval  Space  Surveillance  System 

Dahlgren,  VA  22448 

Library  of  Congress 
Washington,  DC  20540 

ATTN:  Gift  and  Exchange  Division  (4) 

Defense  Printing  Service 
Washington  Navy  Yard 
Washington,  DC  20374 

National  Oceanic  & Atmospheric  Administration 
National  Ocean  Survey 
Rockville,  MD  20850 

National  Aeronautics  & Space  Administration 
Washington,  DC  20546 
ATTN : K.SS-10 

Goddard  Space  Flight  Center 
Geodynamics  Branch 
Greenbelt,  MD  20771 

Aerospace  Corporation 
2350  East  El  Segundo  Boulevard 
El  Segundo,  CA  90245 
ATTN:  Library 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  MD  20910 

Local : 

DK-01 

DK-10  (35) 

DK-ll  (Groeger)  (20) 

DK-12  (Manrique)  (20) 


DK-50 

DK-70 

DX-21 

DX-222 

DX-40 

DX-43  (Elliott) 


(5) 
(3) 
(2) 

(6) 


